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Abstract 

We resolve a question of Bramson and Griffeath by showing that the expected capture 
time of four independent Brownian predators pursuing one Brownian prey on a line is finite. 
Our main tool is an eigenvalue estimate for a particular spherical domain, which we obtain 
by a coning construction and domain perturbation. 

In this paper, we examine the expected capture time of a single Brownian prey pursued by 
n independent Brownian predators. All motion is restricted to a line. Bramson and Griffeath 
|BG| first considered this problem, and estimated the capture time in various circumstances. In 
particular, they showed that if at time t — there are predators on both sides of the prey then 
the expected capture time is finite. For this reason, we will assume that the initial position of the 
prey is iro(O) = 1 and the initial positions of the predators are xi(0) = ■ ■ ■ = x n (Q) = 0. In this 
case, [BG showed that the expected capture time is infinite for n — 1,2, 3, and conjectured that 
it is finite for n > 4 (as indicated by simulations). Li and Shao |LS| showed that the expected 
capture time is finite for n > 5. Using a similar method, we resolve the remaining case by showing 
that the expected capture time is finite for n = 4. 

Theorem 1 Let Xo(t) be pursued by x±(t), . . . , x n {t), where xq, ■ ■ ■ , x n are independent, standard 
Brownian motions on a line, xo(0) = 1, and xi(0) = • • ■ x n (0) = 0. Then the expected capture 
time is finite if and only if n > 4. 

One key difference between Li and Shao's method in |LS| and our method here is that they 
consider a difference process, while we do a geometric splitting in spherical polar coordinates to 
reduce the dimension of the problem. See Section[21for more details on our dimension reduction. 

One can reformulate the capture problem described above in terms of exit of a Brownian 
particle in M. n+1 from a specific cone. We denote the position of the prey at time t as xo(t) and 
the position of the jth predator at time t as Xj(t). By our choice of initial conditions, the initial 
position of the Brownian particle x(t) = (xo(t), x\(t), . . . , x n (t)) is a;(0) = (1,0,..., 0). The event 
of capture is then equivalent to the Brownian particle x(t) first leaving the cone 

Cn+l ■= {{x ,Xi, . . .,X n ) | X > Xj,j = 1, . . . ,n}, 

with a;(0) = (1,0, .. . ,0), and so we must estimate the expected first exit time of a Brownian 
particle from the cone C„+i, with the starting position (I, 0, . . . , 0). 

DcBlassie |DB| developed the theory of estimating exit times for Brownian motion from cones 
in Euclidean space. Let C — {(r, 8) \ r > 0, 9 6 D C §"} be the cone over a domain D C Also 
let t x be the exit time from C of a Brownian particle with starting position x, and let ¥(t x > t) 
be the probability that t x > t. DcBlassie showed ¥(t x > t) ~ c(a;)i~ a< -"- ) , where 

1/2 



2a(n) 



Xx(D) 
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Here \i(D) is the first Dirichlet eigenvalue of D. In the particular case we are interested in, 
the expected exit time of a Brownian particle from C n +x is finite if and only if a(n) > 1, which 
reduces to 

Ai(X>„) > 2n + 2, (2) 

where V n = C n +x nS n . 

Our method for proving Theorem ^ is to estimate the first eigenvalue of T> n using the mono- 
tonicity property of eigenvalues, a coning construction, and domain perturbation. The rest of the 
paper proceeds as follows. In Section^we discuss the geometry of C n+ i, T> n , and related regions. 
Section [21 contains the separation of variables background one needs to estimate the expected 
capture time. For the reader's convenience, we also include the proof that the expected capture 
time is infinite for n = I, 2, 3 predators in this section. We prove the n — 4 eigenvalue estimate in 
Section^ Finally, we describe a numerical computation of \x(T>^) and a lower bound for Xi(T>4) 
in Section E| 

Acknowledgement: Davar Khoshncvisan told us about this problem. Wc thank him and Pedro Mcndcz for lending 
their ears and expertise in this project. A. T. thanks Frank Stcngcr for suggesting the algorithm of Scctionl4l 

1 Geometry of the cone C 

The cone C„+i and its spherical angle T> n have much symmetry. First observe that C n +x contains 
the line C spanned by (1, 1, . . . , 1); this is the line where all the inequalities x$ > xj, j — 1, . . . , n 
are equalities. Thus we can split C„+i as a sum of a line and lower dimensional cone 

C„+i = C® V„, 

where V n = C C n +i- Notice that Vj := eg — Bj, j = 1, . . . , n is orthogonal to (1,1,..., 1), so 
Vx , ■ ■ ■ , V n provide a basis of C x . It is convenient to define 

T n -x :=V n nS n -\ 

where § n_1 is the unit sphere in C = span{V}}. The domain T> n is a double cone over T n -x- 
More precisely, let N be one of the intersection points in S™ n C (there are two such points), and 
let (r, 6) be polar coordinates in S™, centered at N. Then 

-D n = {{r,6) | 0eT n _i,O<r<7r}, 

where we have identified S™ _1 with the unit sphere in the tangent space T^W 1 . In Figure ^ we 
show T~x and T>2 ■ 

In later sections, we will use a generalization of this type of spherical cone. In general, let f2 
be a domain in the equatorial § n_1 of E> n , and let ro £ (0, it}. Then we define the truncated cone 

TC{n,r ) := {(r,6) \0£n,O<r< r Q }. 

We abbreviate TC (Q,tt) = TC(fi). In this notation, V n = TC{T n -x) = TC(T n -x,n). 
The domain T n _x has symmetry. If we let 

C n+1 = {( X 0,%1, ■ ■■Xn) I Xj > X k ,j ^ k}, 

then we see C n+ x = C° +1 and C° n+l are pairwise congruent. Thus T n ~\ is a face of the regular 
(ri + l)-hedral tesselation of the standard one obtains by connecting the vertices of a regular 
(n + l)-simplcx with great circle arcs. In particular, one can compute the diameter of T n -x as 
the distance from a vertex to the center of the opposite face, which is 



S(n — 1) = arccos(— 
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Figure 1: This figure shows 71 and T>%. 



Moreover, the spherical angle of T n -\ at a vertex is T n -2- So we can construct a succession of 
comparison domains for T\ , T2, . . . starting with T% and using the coning construction described 
above. To this end, we let 

t 1 :=[0,y]=T 1 , f n :=TC(f n ^5{n)). 

By induction, T n <zT n , and so \i{T n ) > Xi{T n ). 

2 Separating variables 

We discuss two types of separation of variables in this section. The first type is the separation of 
variables used in |UB| to estimate expected exit times of Brownian motion from Euclidean cones, 
and the second is the separation of variables one performs to estimate eigenvalues of a spherical 
domain with a conical structure. 

One needs the domains in question to have a certain amount of regularity in order to separate 
variables. In particular, these domains are piecewise C , satisfy an exterior cone condition, and 
are well approximated (in terms of Hausdorff distance) by finite volume, C 1 domains. See the 
introduction to DB for a precise statement. It is straight- forward to verify that all of the domains 
we consider satisfy DeBlassie's hypotheses. We will refer to these domains as nice. 

2.1 DeBlassie's separations of variables 

We first review DeBlassie's DB argument. Consider the cone C S W l+1 over a domain flcS": 

C = {r8\ r > 0,6 E D C S"}. 

Let t x be the time it takes for a Brownian particle to exit C, with starting position x, and let 
u(x, t) = P(r x > t) be the probability that t x > t. Then u satisfies the heat equation 

u t = jAu, (x,t) € C x [0,oo); 

u(x, 0) = 1, x € C; 

u(x, t) = 0, (x, t) e dC x (0, 00). 

In polar coordinates (r, 6, t), the PDE becomes 

TL I 

2u t — u rr + —u r + —A§ n u. 
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Moreover, the solution scales as u(r,9,t) = u(f3r,9, /3 2 t), so we can separate variables and look 
for a solution of the form u — R(£)U(9), where £ = r 2 /2t. Then we get a positive separation 
constant Xj(D), and 

, ^ _ A §n U _ A^R + (4£ 2 + 2n£)R 
- ~U~ R ■ 

Letting R = £°p(-£), this ODE becomes 

n + 1 

i'P + ( a H 2 ^ ~ ap = °' 

which has solutions of the form = /(a/2, a + (n + l)/2, £), where / is the confluent hyper- 
geometric function (see e.g., L ). Using / and restricting to K C D compact and T > 0, one 
obtains (see |DBp an expansion for it of the form 

where t/j is the jth Dirichlet eigenfunction of A§n on I? and 



2aj(n) 



1 \ 2 
U — 1 



1/2 

n. — 1 



This sum converges uniformly in K x [T, oo). The leading term in the expansion is u(x,t) ~ 
i3it r i(a;/|a;|)(|a:|/(2t))~ ai (™\ which yields the conclusion of inequality J2J). 

2.2 Separating variables on the sphere 

Next we separate variables to relate the eigenvalues X\(T> n ) and Ax(^-x)- First recall that we 
can write the Laplacian for S" in polar coordinates as 

Agnu = u rr + (n — 1) cotrw r + esc 2 rAgu, (3) 

where Ag is the Laplacian on the equatorial S™ . This lemma gives us a dimension reduction. 

Lemma 2 Let ft be a nice domain in an equatorial S" _1 with first eigenvalue A = Ai(51) , and 
let T> = TC(fl) be the double cone over tt. Then the first Dirichlet eigenvalues of T> and are 
related by 

Ax (2?) = Ax (0) - ^ + ^iH^ + Vfi). (4) 

In particular, X\(T>) > 2n + 2 whenever Xi(Q) > 2n. 

Proof: Set w(r, (9) = R(r)T(9), where i?(0) = = #(tt) and T(0) = for 9 € 517. Then u is 
an eigenfunction on I? with eigenvalue /i precisely when 

TR + (n— 1) cot r TR + esc 2 r i?A e T = -fj,TR. 

Separating variables with a positive separation constant A yields 

sin 2 rR + (n — 1) sin r cos rR + p sin 2 ri? A#T 

= A = . 

R T 

Choosing T to be the first eigenfunction of £1, we obtain the ODE 

sin 2 rR + (n - 1) sin r cos ri? + (p sin 2 r - A)i? = 0, (5) 
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which has regular singular points at r = 0, 7T. If we try a solution of the form R = sin" 1 r, for 
some power to, we find 

= m(m — 1) sin m r cos 2 r — to sin m+2 r + m(n — 1) sin m r cos 2 r + (jjt sin 2 r — A) sin™ r 
= sin m r[(m 2 + m(n — 2) — A) cos 2 r + (/i — m — A) sin 2 r] . 

Both coefficients must vanish, so we have /i = m + A and A = to 2 + m(n — 2). Solving for to, we 
find 

2-n /(2-n) 2 
m = — + V 4 +A - (6) 

Next, observe that if Ai(f2) = 2n then Ai(2?) = 2n + 2. Finally, the formula for Ai(2?) is monotone 
increasing in Ai(Q), and so Ai(2?) > 2n + 2 whenever Ai(f2) > 2n. □ 

Remark 1 ^4 second solution to equation has the form sin™ r cos r, where to is again given 
by equation 0) but 

H = A + 3m + n. 

TTizs eigenfunction vanishes on {7r/2} x 3 so it corresponds to a higher eigenvalue. 

At this point, we can prove that the expected capture time for n = 1, 2, 3 predators is infinite. 
To prove the expected capture time is infinite, by inequality (J2J we need to show Ai(2?„) < 2n + 2, 
or, equivalently, that Ai(7^_i) < 2ri. In the case of n = 1, we have T)\ = [— 37r/4, 7r/4], and so 
Ai(2?i) = 1 < 4. In the case n = 2, we have T\ = [0, 2n/3], and so Ai(Ti) = 9/4 < 4. We cannot 
compute Ai(72) so easily, but we can find a test function to show that Ai(72) < 6. We shall find 
Ai(72) numerically in Section 0] Recall the Rayleigh characterization of the first eigenvalue of a 
domain Q: 

'SaW 



Ai(fi) = inf 



r2 



To show that Ai(T 2 ) < 6, it suffices to find f € #£(7^) so that / Ta |d/ | 2 / j % f$ < 6. Let 

fo(x) = sin(dist(x,<9T 2 )), 

and observe that |d/o| 2 = 1 — /o off the set of focal points of 972, which is a set of measure zero. 
A computation shows 

Ai(T 2 ) < = < 6. 

Jr 2 ^0 tt - V3 

One can generalize the eigenvalue relationship (0J to spherical cones of the form TC(fl,ro), 
for < ro < 7r, using the hypergeometric function 

, . _ a(a + l) z 2 a(a + l)(a + 2) z 3 

/(a, (J, j,z)- +--+ ^ - - ^ + ^ + + 2) ^ - — + . . . () 

Later we will use the following lemma to relate the eigenvalues of T n —x and T n , which we will 
then estimate to bound the asymptotic decay rates of probabalistic exit times. 

Lemma 3 Let be a nice domain in a equatorial S™" 1 C §" with eigenvalue A = Ai(f2). Then 
the first Dirichlet eigenvalues of T2 and TC(il, r ) are related by 

Ai(TC(Q,r )) =/i = fi(n,X,r ), (8) 

where fi is the first Dirichlet eigenvalue of the ODE Jj5J) on [0, ro] . Ifro>n/2 then fi is the unique 
zero of f(ai, f5\, 71, (1/2)(1 — cosro)) in (to + A, 3to + A + n), where to is defined in equation Jj^) 
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and 

1 + J(n- 2) 2 + 4A ± J(n - l) 2 + 4^ 
ati,Pi = - - 



2 + J(n - 2) 2 + 4A 
71 = 2 ■ 

-(f r o < </ien /j. «s i/ie unique zero of f (012, /?2)72 , (1/2) (1 — cosro)) in (0,n), where 

n - 1 ± J(n - l) 2 + ifj, n 
a 2 ,P2 = 2 ' 72 = 2' 

with f defined by equation Q). 

Proof: We separate variables and look for a solution of the form R(r) = sin™ ru(r), with 
u(r) 7^ on [0,r ), but u(r ) = 0. Then equation JSJ) becomes 

= sin m+2 ril + (2m + n — 1) sin m+1 r cos r u + [m(m + rt — 2) cos 2 r + (/1 — to) sin 2 r — A] sin" 1 r w. 
Now let u(r) = y(x), where x = (1/2) (1 — cosr), which transforms the ODE above into 

x{\ — x)y" + (m + — n — (2m + n)x)y' — (A + m — fi)y = 0. 
The solution to this ODE is the hypergeometric function y(x) — f(a,[3,-f,x), with 



2to + n - 1 ± J (2m + n- l) 2 - 4A - 4m + <kt 
a,/3 = 3 

2m + n 

7 = —2~- 

The lemma follows from taking R(r; n, A, ro) = sin m r /(a, /3, 7, (1/2) (1 — cosr)), where we choose 
\i so that i?(ro; n, A, ro) = 0. □ 
One can use this lemma to compute Ai(7^,_i) iteratively. In this case we start with Ai(Ti) = 
9/4 and apply equation (jHJl. For the reader's reference, we include a table of \\{T n -\) and the 
corresponding lower bound for a(n) (see equation l|T|l). In this table we use equation (0J and the 
fact that V n c TC(f n ^ 1 ) so Ai(X>„) > \ 1 (TC(f n ^ 1 )). 

n Xi(T n -i) lower bound for a(n) 

2 2~25 .75000000 

3 5.00463581 .89614957 

4 7.884040724 .99030540 

5 10.77018488 1.05417466 

6 13.6203196 1.09882819 

From this table, one can see that X1CT3) is slightly less than 8, so T3 comes close to, but is not 
quite, an effective comparison domain for proving Theorem^ 



3 The eigenvalue estimate 

In this section, we complete the proof of Theorem ^ First observe that it suffices to show 
Ai(7ji) > 8. In fact, it suffices to find an effective lower bound for Ai(7^) using the following 
scheme. We define A cr by 

8 = /i(3, A cr , <5(3)), A cr ~ 5.101267527. 
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If we construct a domain Q2 C § 2 such that T2 C Q2 and Ai(0a) > A cr , then by domain mono- 
tonicity and equation JSJl, 

Ai(T 3 ) > A^TC^,^)))) > Ai(TC(e a ,«(3))) > 8, 

completing the proof of Theorem ^ As the last step in our proof, we construct a domain Q2 C § 2 
as a perturbation of 72, such that T 2 C & and \\[^%) — 5.102 > X cr . 

Figure El shows sketches of the domains T2, Q2 and T2, after stereographic projection to the 
plane. 
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Figure 2: The inner dashed curve is T2, solid curve is Q2, and the outer dashed curve is Ty- We 
generated this figure using the computer program MAPLE. 



Previous work of Rayleigh [R] and Polya-Szego |PS| motivates us to consider this type of 
domain perturbation. They studied the eigenvalue of a planar domain which has the form 
{(r, 9) | < r < c + ef(9)} in polar coordinates, for some small e > 0, and gave an expres- 
sion for the eigenvalue Ai in terms of ef. In our case, we fix A = Ai and find a domain Q2 with A 
as its first eigenvalue. 

Suppose the functions R(r) and 0(9) satisfy 

e" + A6 = 0, 0<6><^; 

9(0) = 0, 
9(2tt/3) = 0, 

sin 2 r R + sin r cos rR+(fi sin 2 r — X)R = 0, < r < tt; 

R(0) = 0. 

Then u(r, 9) = R(r)Q(9) is the first eigcnfunction of A§2 on its nodal domain G2- By construction, 
A§2it + fiu = 0, so u is an eigcnfunction. Also, u does not change sign on its nodal domain, so it 
must be the first eigenf unction. In polar coordinates, u(r, 0) = = u(r, 2ir/3). Let m = \f\ and 
set R(r) = sin™ ru(r), so that equation (J5J becomes 

sin 2 (r) il + (1 + 2m) sin(r) cos(r) u + (/1 — m — \)u = 0. 
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Next we take A = 9/ 2 /4, corresponding to the I th mode of the interval [0,2tt/3], and write the 
solution in terms of the hypergeometric function: 

ui{r) = /(3Z/2 + .5 ± v/l/4 + /i, 1 + 31/2, (1 - cosr)/2). 

Finally, we take the fi = 5.102 superposition of the 1 — 1,3 modes to define 

sin 3/2 (r) sin(30/2)/f(r,0) := (sinr) 3/2 ui(r) sin(3(9/2) - .0003(sinr) 9/2 w 3 (r) sin(9(9/4) 

and let G2 be the nodal domain of sin 3/ ' 2 (r) sin(36'/2)iJ. By construction, Ai(^) = 5.102 > X cr . 

The domains T%, Qi and Ti all have three boundary curves, and two of these boundary curves 
for each domain lie along the great circle arcs 8 — 0,27r/3. For convenience, we convert to 
a planar domain using stereographic projection, with the south pole corresponding to r = 
in polar coordinates on § 2 . The relationship between polar coordinates (r, 8) on S 2 and polar 
coordinates (p,6) on the plane is p = tan(r/2). Under stereographic projection, the great circle 
arcs 8 = 0, 27r/3 correspond to rays emanating from the origin at angles 0, 27r/3. The third 
boundary curve of T2 is given by the circular arc 

( ^) cose __L)2 + w) _y| )2 = | 

which we can rewrite as 

a . n . y/2 cos(6» - tt/3) + ^2cos 2 (6»-7r/3)+4 
P = (3(8) = 1 ■ 

Thus showing T2 C Q2 is equivalent to showing H (r, 8) > along the arc given by (p = (3(8), 8) 
for < 8 < 2tt/3. 

One can see from Figure |21 that H > along this third boundary component of 7~2 ■ One 
strategy for a rigorous proof that H > is the following. We first evaluate at a point do on the 
curve, checking H > at (p — (3(8q), 8q), and bound the derivative Hg on an interval containing 
8q. Our bound Hg > M gives us a lower bound H > H(8q) — M(6 — 6q). Thus H > on a possibly 
smaller neighborhood of 8q . We then repeat this procedure with each endpoint of this (smaller) 
interval. One can simplify the computations by observing T 2 and Q2 are symmetric about the ray 
8 = n/3. Then it suffices to evaluate H and its derivative at 8 = 0, 1/2, 2/3, 27r/9, 7r/3. □. 

Given our lower bound Ai^) > Ai(<?2) = 5.102, one can use equation © to show AiCTjj) > 
8.00087815. This in turn gives Ai(D 4 ) > 10.001024501 and a (4) > 1.00007318. In the same way, 
the lower bound Ai(T 2 ) > 5.102 gives a(3) > .90671950. 



4 A numerical computation of the first eigenvalue 

In this section we describe a numerical computation approximating the eigenvalue \\(T2), which 
relies on Stenger's sinc-Galerkin scheme [St]. We will show \\(T2) » 5.159..., and so a(3) « 

.9128 .... Also, using equation ©, we have Xi(T 3 ) ~ 8.0691 and so o(4) ~ 1.0057. 

Given h > and a positive integer fc, one defines the fcth cardinal sine function of stepsize h 
as S(h, k)(hk) = 1 and if z 7^ hk, 

Following StengerjSt], for z = a; + iy we shall let the basis functions be 4>jk(z) = ay(#) x (3k(y), 
where 

n 

aj(g) = 5(j, fe) o ln( ^ ), a n+ i(x) = sin 2 (x) - sin 2 (xi)ai(x) (10) 
2 _ x i=_„ 
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for 7 — —n . . . ,n with sine points xi — — rrr and 

2(1 + e hl ) 

n 

(3 k (x) = S(k, h) o logsinhy, (3 n +i(y) = sech(y) - ^ sech(yi)0i(y), (11) 

l=—n 

for k = —n, . . . ,n with sine points yi = sinh _1 (e' 1 '). We let x n+ \ = tt/2 and 2/n+i = 0. The 
dimension of the space of functions 4>jk ism = (2n+2) 2 . Also a n+ i(7r/2) = 1 and /3 rl +i(0) = 1, but 
both are zero at the other sine points because their interior sine approximations are subtracted. 
The functions <j)jk = oij x [3k form a basis of sine functions when we are working on a rectangular 
strip. 

The most important property of sine functions is that they are well-suited to approximating 
integrals in a strip. In particular, one can increase accuracy in a numerical computation by 
increasing the number of evaluation points, without recomputing grids (as in finite elements). 
In order to take advantage of this simplicity, we will conformally map a subdomain of T~2 to a 
half-infinite strip. 

Recall that T2 is an equilateral triangle on § 2 with all its interior angles equal to 27r/3. Let 
Ti , T2 , T3 be the vertices of T2 , let Sj be the midpoint of the side opposite Tj , and let F be the 
center of mass of T2 . Observe that T2 is invariant under reflection through the lines FTj , which 
divide 7~2 into six congruent subtriangles (see Figure |3J). Thus the first eigenfunction is invariant 





F 






£2 









S3 



Figure 3: The three marked points are F, T\ and S^. 



under these reflections, and we can recover it by restricting to the smaller triangle fi, which has 
vertices F,T l7 S 3 . The triangle SI has a right angle at the vertex £3, and an angle of n/3 at the 
vertices F and TV The first Dirichlct eigenfunction of T 2 , restricted to fl, will have Dirichlet data 
on the edge T1S3 and Neumann data on the edges FTi and FS3. 

We first transform Q to a planar domain (which we again denote as f2) using stereographic 
project, sending S3 to 0. This transformation sends F to ^(y/E—y/2) and T\ to |(v6 — v2). Now 
SI is bounded by the two straight line segments joining to — \[2) and to | (V6 — \/2) 

and the circular arc joining ^{V® — V%) to ^(V& — \/2) which makes an angle of ir/3 with the 
axes. Next we find a conformal transformation / which takes Q to the half-infinite strip 

L> = {zeC|0<5Rz< tt/2, < 3z} 

such that /(0) = 0, /(| (a/6 - V2)) = tt/2, and f(±(V6- \/2)) = 00. Under this transformation, 
the first eigenfunction u satisfies 

0, lizeD, 

0, if !Jie z — and 3m z > 0, 

7T 7T 

0, if < 5Rc z < — and 3m 2 = 0, or if 3?e z — — and 3m z > 0; 



A « + = 
u = 
du 
dn 
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where A* = /*A = A z is the pulled back Laplacian. 

There are formulae for conformally mapping domains bounded by finitely many circular arcs 
[N] generalizing the Schwarz-Christoffel formula. Following Forsythe [0, the Schwarz triangle 
mapping z G D or from sin 2 z G {Qz > 0} of the upper halfplane to w G Q is given by 



cos 2 z 



(w 4 + 2^3? 



If 



(w 4 -2x/3w 2 -l) 3 ' 



The group generated by reflections of fl along its edges in fact tiles the sphere. The inverse 
function w h- ► z 2 extends to a single sheeted cover and is invariant under the symmetry group, 
thus is a rational function of the plane. Thus we may compute /. Writing g(z) = cos 2 / 3 z, 



/(*) 



1-9 



V / 3(l+. 9 ) + 2 v /l + .g + .g 2 
Pulling back under w = f(z), the conformal weight takes the form 



4 
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dz 
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\-g)+2y/l + g 


+ g 2 




\g\( 


V3(l + g)- 


^2^/1 + g + g 2 


+ |i 


~g\) 2 



The branch cuts for the square and cube roots may be taken above the negative real axis. Thus 
17(D) lies in the fourth quadrant so the denominator in / is nonvanishing. 
Let G(z, z') denote the Green's function for the problem on D 



A* u = g, if z G -D, 

u = 0, if Sfte z — and 3mz > 0, 



du 
On 



0, if < Sftez < 7r/2 and 3mz = or if 5Rez = tt/2 and 3m z > 0. 



The Green's function may be found by the method of images. Thus, denoting w = sin z = 
x + iy, w* — —x + yi and uj — sin^, we get w* = (w)*. The Green's function is 

G(z; () = — (\n \w — uj\ — In \w* — u\ + In \w — uj\ — In \w* — ui\j . 

Pulling back by /, the eigenvalue problem for the triangle may be restated as finding an eigenvalue 
for the integral operator 



u{z) 



G(z,z')\df(z')\ 2 u{z')dz' 
(l + l/(^)| 2 ) 2 



=: Au{z) 



(12) 



The key point is to notice that the operator has logarithmic and algebraic singularities at the 



points z' = 0,7r/2 and 



Thus we need to handle these singularities. The solution is zero 



along the imaginary axis, but free along the other two sides, which may be extended to functions 
to the plane which have odd reflection symmetry along the imaginary axis and even reflection 
symmetry along the other sides. We shall approximate u(z) in an m-dimensional space X m with 
the same symmetries. Also, noticing that the eigenfunction on T2 behaves like 1 — dist 2 (,z,T2) at 
the vertex F, we actually have u G Lip(_D) and u decays algebraically at 00. Choosing a basis 
{0i, . . . , cf> rn } of X m , we shall compute the matrix of the transformation Aik = VeA(f>k, whose 
largest eigenvalue \x m — > 1/A as m -> 00 and which is an upper bound A < 1/ /z m [ot] . The 
integral operator shall be computed numerically via sine quadrature, which can handle such mild 
singularities. 
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We make the approximation u(z) « Vu{z) = W <j)jk ( z ) where the sums are over j = —n . . . , n+ 
1 and k = 0, . . . ,n+l. Here 4>jk — otj xftk, where ct,- and [3k are defined in equations i|l(J[l and (fill) , 
respectively. We define V = VjkU — u(xj + iyu)- Thus the approximation Vu is a collocation, 
as it equals u at the sine points. Thus, the matrix is approximated by 

Ajk, pq = / G{xj,y k ,£,r))(j>p q {£,r})^{£,ri)d£dr) 

J D 

w v^w K G(xj, y k , x t , y^pqjx,, y K ) *(x t ,y K ) 

where '!'(£, 77) = 4(1 + \f{£ + ir))\ 2 )~ 2 \df(( y +iji)\ 2 . The approximating sum is carried over 4m 
terms, corresponding to sine quadratures in the four regions bounded by singularities (e.g. in 
case < Xj < ir/2 and < yk): Di — {£ + if] : < £ < Xj, < f] < yk}, Du — {£, + irj : Xj < 
£ < 7r/2, < T] < y k }, Dm = {£ + iri :0 < £ < Xj, y k < y} and D IV = {S,+ir]: Xj < £, y k < v} 
and where v L and w K are corresponding weights and x L and y K are corresponding sine points. 

The computation shows that XxiTz) » 5.159... which is above the critical value. This is 
computational evidence that the exit time has finite expectation. We present a table of a few 
the computed eigenvalues for approximations in spaces of given dimension, coming from the 
Sinc-Galerkin collocation scheme described. 

Dimension Eigenvalue Estimate 
16 5.948293885960918 
100 5.262319373675790 
1024 5.153693139833067 
2116 5.158585939808193 
2304 5.158832705984016 
2500 5.158849530710926 
2704 5.158968860560663 
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